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COMPABISON OF THE ' ENERGY; METHOD ■ VTTH THE ACCELEBOMETEB METHOD 

■ . \ 

OF COMPUTIN'* DBAG COEFFICIENTS FBOM FLIGHT DATA 

. 1 ■ * 

By Thomas L. Keller and Hohert F; Keuper 

■ 1 1 •» 

, . / summary \ 1 ' 

A comparison was made between the energy and the accel- 
erometer methods of computing drag coefficients from flight 
tests of . t he F— ' 51B and the F—39N airplanes. The energy 
method was found to he unreliable with the present types of 
instrumentation except under the practically unattainable 
conditions of very steady flight at constant indicated air- 
speed within a vertically stationary air mass. 


INTBODUCT I ON 


In the pabt, many measurements of the drag of airplanes 
in flight have involved the use of some variation of the so- 
called "energy method" (described in the appendix) for the 
reduction of the data. The results obtained by this method 
are frequently very erratic, occasionally even giving -nega- 
tive drag- coefficients. 

Since the results from the "accelerometer method" 
(described in' the. appendix) appear to be’ consistent with 
computed and wind-tunnel data, the flight— test results re- 
ported In ref erence .1 . of f er -an excellent opportunity for 
comparison of the two methods to determine the oause of the 
discrepancies. This report presents -such a comparison and 
includes supplemental data obtained from two special divos 
performed at constant indicated airspeed- in the F-39— N-l 
airplane. 


2 


. ' HACA CB No. -SH31 


■ -DZSCBIPTION OF EQUIPMENT . 

Airplanes 

The airplanes UBed for the tests were the North American 
P-51B~ih-NA and the Bell F— 39N-1. Complete descriptions of - 
these airplanes are given in references 1 and 2, 


Instrumentation 

Standard NACA photographically recording instruments, 
installed in both airplanes, were used to obtain airspeed, 
pressure altitude, normal and longitudinal accelerations, 
and (for the P— 39N airplane only) angle of attack as func- 
tions of time. Tallies of. free-air temperature .at different 
altitudes were obtained during ascent from the pilot Vs indi- 
cating instrument, and were corrected for temperature rise 
due to impact of the air. stream. 

. The airspeed— moasuring systems utilised., as sources of 
pitot— static and total pressures, freely swiveling pitot— 
statio heads mounted on booms located beneath the wing tips. 
These booms extended approximately 0.8 of the local wing 
chord ahead of the leading edge. Each pitot— static head 
consisted of two static— pressure tubes and one total— .pres- 
sure tube, which permitted the use of independent sources of 
static pressure fcr both the airspeed and the altitude 
recorders, 

'The pressure lines from the- pit ot— static head to the 
recording instruments were. made ns short as possible in 
order to minimize lag, and the lines tp the recording air- 
speed meter were balanced to give equal flow rates in the 
static and total— pressure tubes during rapid changes in - 
pressure,' Ground tests of nockr-ups of the airspeed and 
the altitude pressure lines indicatod that .the lag in the . 
systems at the maximum rates of descent caused an error in 
the rocordod altitude of about 20 feet at 15,000 feet 
(negligible) for both airplanes. 

A correction for the position error of the pitot— static 
tube was determined from flights at various airspeeds past 
a fixed visual reference point of known constant pressure 
altitude. It was assumed that the measurements of total 
pressure were correct and that the variation of recorded 
altitude with airspeed (obtained by a special sensitive 
altimeter) at the constant pressure altitude resulted from 
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the position of the static tube only. Since the maximum 
error in altitude^ ss determined by this calibration, was 
loss than the accuracy of - the ■ recording-' dltimet er used in the 
dive tests, no attempt was made to -correct the altimeter read- 
ings for position error. 

The accuracy of the swiveling pitot— static head has been 
investigated at Mach numbers up to 0.80 in the. 16— foot wind- 
tunnel at the Ames Aeronautical ■ Laboratory . The results Bhow 
that the effectB of - compressibility not. included in the regu- 
lar airspeed calibration ar £ negligible over the flight Mach 
number range investigated. 

The correst Indicated airspeed was computed by. use of 
the standard formula 

V t = 170 3 |( 5=2 ♦ i y-°° 6 _ ! ] w * 

, L ' P 0 ' 

whore 

correct indicated airspeed, miles per hour. 

H. free— stream t otal ’ pre s sure 

i 

p free— stream static pressure 

p Q standard atmospheric pressure at sea level 


METHOD OF ANALYSIS 


If accurate drag data are available for an airplane in 
flight it is possible to ahalyse, indirectly, the errors 
involved in the use of the energy method by working backward 
from the drag data, through the energy method equations, to 
obtain values of airspeed or altitude. A comparison of these 
values with tie actual values ' recorded during’ the flight 
offers some insight as to the nature of the errors. The fol- 
lowing discussion gives a derivation of the necessary equations 
and a detailed explanation of' the process. 

The equations .for the drag coefficient, "based on the 
energy and the acceleroaeter methods (excluding propeller 
thrust) have been derived in the appendix. The equation from 
the energy method has been shown to be 
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if ( dh/dt 




( 1 ) 


g 


■whore* ■ 

■; ' 

V " weight of a-irpla.ne ,■ pounds 

q free— strerm dynamic pressure, pounds per square foot 
S wing area, square feet 

h true altitude, feet 

« . j 

•» * '* 

true airspeed, feet per second 
t time, seconds 

g acceleration due to gravity, 32.2 feet per second per 

second 

In this discussion the slope of the altitude curve 
dh/dt will he assumed postive when h decreases, and the 
slope of the airspeed curve dVip/dt will he assumed posi- 
tive when Vip increases; hence the minus sign in equation 
(1), Equation (l) .can he transformed as follows! 

dVfj, f dh/dt CpqS^ 

4t = J * 


( 2 ) 


or 


dh / c p 
dt \ W 


C D qS dVj/dt 




(3) 


The . f ollowing expressions, solvable by graphical in-r 
tegration, . can he derived from equations (2) and (3): 

... 7- T m . + g it 

T . T 0 J 0 ' 7 T W ) 6 4 

pt /o^qS dV-r/dt \ 

h o h 0 - / ( -Jd- + _i_ ) t t at 

■Jo V W e ' 
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Thus-, if the drag .date obtained by the accelerometer ■ 
method are approximately correct and the reaorded values . of 
altitude are assumed to be correct, the corresponding - air- 
speed— time curve can be computed by a series of approxi- 
mations and integrations. In the same manner, the corre- 
sponding altitude— t ime curve can be computed if the recorded 
values of airspeed are assumed. to be aorrect, The dif- 
ference ■■■ between • the computed, and recorded airspeed or alti- 
tude- curves-,., thep;. is th®-- amount of error which' would pro- 
duce the corresponding, error in ' the drag-coefficient curve, 

' TESTS, RESULTS,. AVD DISC.USSIOV 

Unrestricted Plight Conditions 


Figures 1 and 2 show examples of airspeed and altitude 
calculated by the foregoing method from the data of a 
propeller— off dive of the B-51B airplane (dive Vo. 2, 
reference 1). Figures 3 and 4 show drag curves for the same 
flight, calculated by the energy and the accelerometer method 
from which it 1 b apparent that the energy method yields 
extremely erratic results. 

The errors in the altltudo curve (fig, 1) or in the 
airspeed curve (fig, 2) necessary to produce the variation 
of drag coefficient shown in figures 3 and 4 are very small, 
and it is apparent that a small error in either airspeed or 
altitude could cause sufficient change in slope to produce 
a very large error in the computed drag coefficient. In 
fact, evon small errors in fairing these curves could lead 
to appreciable errors in drag coefficient. Furthermore, it 
should be noted that the curves of figures 1 and 2 represent 
the outside limits of the error;- that is, the error in each 
quantity was computed on the basis of no error in the other 
quantity. Actually, the error probably is divided between 
airspeed and altitude and thus is leas than that shown for 
either. 

The dives reported in reference 2 did not follow any 
definite flight technique; that is, airspeed and altitude 
were allowed to vary linr estr ained, the sole objective of the 
pilot being to attain as high Mach numbers as possible. The 
question has been raised, however, as to whether dives at 
constant indicated airspeed would yield data amenable to 
the energy method. It was believed that this type of dive 
would yield the smoothest possible variation of true airspeed 
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with t-ime, vi.th the resultant increase in accuracy of the 
-measured slopes dYg/dt and decrease in any unknown errors 
'.due to rapid pressure changes.. 

*■ » , 

RSstricteH flight- Conditions 

In order to ascertain the effect of a definite flight 
technique .on the results of the energy method, therefore, 
two dives at nearly constant indicated airspeed were per- 
formed in the P— 39B-1 airplane. A sensitive longitudinal 
accelerometer was includod in the instrumentation for these 
flights, so that a direct comparison could he made between 
the doubtful quantities involved in tho energy method and 
corresponding quantities of established accuracy used in the 
accelerometer method. 

Time histories of variables measured during the two dives 
are presented in figures 5 and 6. The necessary corrections 
for free— air temperature were applied to the Indicated values 
of airspeed and rate of change of pressure dp/dt to obtain 
the true values of airspeed and rate of change of altitude 
dh/dt used in the computations of the longitudinal foroe 
factor. 

figure 7 shows a comparison of the longitudinal force 
.factors (Ag sin a — Aj cos. a) from the accelerometer method 

from the energy method, ' Theoretically, 

during any one dive, these quantities should be identical. 

The resultant drag coefficient computed from these data is 
presented in figure 8. It will bo observed that during the 
time interval for each dive in which the conditions were 
extremely steady (70 to 145 seconds in dive number 1, and 
55 to 150 seconds in dive number 2) the maximum difference 
in drag coefficient ACp Is approximately 0.0021 and the 
average ACp is very small. The agreement in results is 
considerably improved over that of the unrestricted dive. 

(See figure 3.) It is apparent, however, that the slightest 
variation from st eady. ,c ondi t i on s causes wile variation be- 
tween the computed longitudinal force factors and renders 
the energy method unreliable. 

furthermore, . since tho derivation of the energy method 
does not allow for any. extraneous acquisition of energy, 


dh/dt dV„/dt 
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such aa that duo to vertical air ourrents or changes of 
vert ical air currents,' c onsiderable error -could he intro— 
auced hy the airplane's transit through a rising or de- 
scending alr,-masB.- - * . 

It is beLiaved that; if extremely steady flight condi- 
tions (including constant indicated airspeed) could he main- 
tained hy the pilot, the energy ■ method of computing drag 
coefficient would he of value. 1 In practical f 1 ight— res earch 
work, especially at high Mach .numbers, however; the required 
conditions are seldom achieved and are difficult to recognize 
when they are achieved.. Thus,, i-t is believed that the 
energy method cannot he relied upon for conclusive results. 

CONCLUSIONS 

From comparison with the accelerometer method, the 
energy method of computing drag coefficient from flight data 
is believed to he unreliable with the present instrumentation 
for the following reasons: 

1. The conditions of constant indicated airspeed and 
very steady flight, which appear to he essential for accurate 
application of the energy method, are impractical for ordi- 
nary flight— research work (especially at high Mach numbers) 
and are difficult to achieve and to recognize. 

2. The energy method is correct only for flights in air 
having no vertical velocity, 

j^mes Aeronautical Laboratory 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 

. x It is instructive to note that, since, theoretically, 
the energy method Bhould he exact for flight in still air, it 
is possible that tho altitude and airspeed curves of figures 1 
and 2 represent a type of pressure calibration and that the 
recording systems are actua-lly in error hy various amounts. 

The differences shown in figures 1 and 3, however, only can 
give a qualitative picture of the errbr because the relative 
amounts of the total error attributable to altitude and to air- 
speed are not known and tho vertioal motion of the air mass 
during the dives could not be determined. Further study along 
this line might be of considerable interest. 
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•APPENDIX A 
Symbols 

, i , , 

1 Symbols used throughout .this report. are listed as follows: 
Vi correct indicated airspeed, miles pier h'cur 
E free— stream- total’ pressure . 
p . free— stream static pressure 
p c standard atmospheric pressure at sea level 
W weight of airplane, pounds 

q ' free— stream dynamic pressure, pounds per square foot 
S wing data, square feet 

h . .true altitude of airplane, feet 

. i 

Vj . true airspeed, feet per socqnd 

t * ■ • 

t time, seconds 

g acceleration due to gravity, .32,2 feet per second per 

second 

E energy, foot-pounds 

D drag, pounds 

Ag the ratio of the net aerodynamic force along the air— 

. • plane. Z— axis (post.ive. when direct ed, upward), to the 
weight 'of the /airplane 

r - . . 

■ • * # . ...»: t . ■ . 

Aj. the. ratio of. the net aerodynamic force along the air- 
• . plane. X— axis (positive when directed .forward.) to 

the weight of the airplane 

< 1 ■ 

E the resultant .of A% and Aj . 

■ ’ ■ if,.*.. 

ft flight— path angle from hor 1 z out. al ,■ degrees 

a angle of attack, degrees 


I 
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APPEND IX B 

DERIVATION OF TEE ENERGY AND THE ACCELEROMETER METHODS 
OF COMPUTING DRAG COEFFICIENT FROM FLIGHT DATA 
Derivation of Energy Method 


Tho total energy (excluding that due to propeller thrust) 
of an airplane In flight in still air with respect to a point 
on the earth at sea level may he expressed by the eauatlon 

w(v T ) a 

E = Wh + ■= (Bl) 

2g 

The rate at which energy is expended by n, n airplane in flight 
is equal to the drag multiplied by the velocity, or 

[l ' ' ' 

AH e DT - y. ( B3 ) 

dt \ 1 

Differentiating equation (31) with respect to time gives 

£5 = + II ££t\ ( B3) 

dt / 

(B3) yields 
I? <jVT >| 

6 dt J 

+ 

8 ' 


+ ilxLil) 

XT ' 


dt ^ dt g 

Combining equations (B2) and 


or 


■ D -~ ¥ u? + 
»-* (iWii 

v T i 


and thus 


CD 


JL ( dh / 

qS V Y 


dt 
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Derivation of Accelerometer Method 

The diagram of the accelerations (excluding that due to 
propeller thrust) acting on an airplane in flight 1 b given 
bc-lows 



The resultant acceleration acting on the airplane is R, 
Its components along the airplane axes are Ag and A z , 

The external (aerodynamic) forces which produce Aj 
and A Z are equal to WAg and WAg and act in the same 
diroctions as Aj and A z »- 

Poeitive drag is defined as an oxtornal force acting 
rearward along the flight path. Therefore, since the com- 
ponent of positive WAj acts forward along the flight path 

and the component of positive WAg sets rearward, for posi- 
tive, values of a, the drag equation "becomes 

D = WA Z sin a — WA X cos a 

and 

°D = ( A Z sln c - A X coa o) 


HACA CB No. 5H31 


11 


BBP3BENCES 


1, Nisson, James M. , fladebsrg, Burnett, L. , and. Hamilton 

William T. ; Correlation of the Brag Characteristics 
of a P— BIB Airplane Obtained from High-Speed Wind— 
gunnel and Plight Tests. VACA ACH Ho, 4K02, 1944. 

2. Gasich, Welko 3., and Clousing, Lawrence A.} Plight 

Investigation of the Variation of Drag Coefficient 
with Mach Number for the Bell F— 39N— 1 Airplane. 

NACA ACH Ho. 5D04, 1945. 





NACA CB No. 5H31 


rig. a 



14 18 88 26 30 34 38 48 46 60 54 

Time, «»o 


figure 2.- Time history of airspeed and acceleration along the flight path showing comparison of 
values recorded and calculated by the energy method. North American P-51B-1 airplane 
with propeller removed. 


Aoaeleromet er method 
Energy- method 



Figure 3.- Time history of drag coefficient ae calculated by the 

acceletoaeter and energy methods . North American P-51B-1 
airplane with propeller removed. 


Figure 4.- Variation of drag coefficient with J£ach number as 
calculated by the accelerometer and energy methods. 
North American P-51B-1 airplane with propeller removed. 
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True airspeed, Vt, fpe Static pressure 



Indicated airspeed, Vj., mph 


True airspeed, Vt, fpB 



figure 6,- Time history of dive No. 2 In which the pilot attempted to hold constant Indicated airspeed Sell P-39N-1 airplane, engine 
throttled, propeller in high pitch. 
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figure 7.- Variation of longitudinal force factors from - accelerometer and energy methodB for two dives at approximately constant indicated 

airspeed. Bell P-39N-1 airplane with propeller in high pitch. -j 
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Reoultant drag coefficient 



Figure 8.- Variation of the resultant drag coefficient calculated by the energy and accelerometer methods for two dives at approximately ? 

constant indicated airspeed. Bell P-39H-1 airplane with propeller in high pitch. a 
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